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A mathematical model of heat transfer in a firestopped vertically 
oriented cable bundle is derived to assist in planning fire test exper- 
iments and to enumerate how changes in geometric and thermophys- 
ical properties affect the temperature rise in the cables when subjected 
to standard furnace fire tests. The analysis indicates that the primary 
heat transfer mode to the cable array is from the flow of hot furnace 
gases up through the void spaces between the individual cables. As 
expected, the most practical and effective way of reducing the heat 
transmission characteristics of a cable bundle is by tightly packing 
the firestop, which reduces the void space between cables and provides 
heat sinking to the cooler environs. 

I. INTRODUCTION 

In this paper a mathematical model of heat transfer in the cable 
bundle of a firestopped vertical cable assembly is developed to assist 
in planning experiments and to evaluate the relative influence of the 
geometrical and thermophysical properties of this portion of the fire- 
stopped configuration. A representation of this complex cable bundle 
geometry in terms of an approximate transient one- dimensional, 
lumped parameter model is obtained through heuristic arguments. 
This is accomplished in a systematic way by first deriving a simplified 
model for a single conductor wire and progressing up in scale, by 
averaging and lumping parameters, to arrive at a heat transfer model 
of a single cable. The heat transfer in a cable bundle is then treated 
using the single-cable model. Cable-to-cable heat transfer is handled 
through boundary conditions at the contact surfaces of the individual 
cables. A set of coupled transient one- dimensional equations results, 
with as many equations as cables in the array. The model is then 
exercised to compute transient temperature distributions within a 
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Fig. 1 — ASTM El 19 standard temperature vs time curve. 

vertically oriented firestopped cable bundle when exposed to the 
ASTM E119* temperature variation shown in Fig. 1. 

Properties of three generic-type cables (terminating, switchboard, 
and power cable) are used to suggest how physical and geometrical 
characteristics of the cables influence the effectiveness of a firestopped 
cable closure. It should be noted that an absolute evaluation of a 
firestop using the approximate heat transfer model developed here is 
beyond the capability of the model. This can only be attempted when 
considering potential nonlinear combustion modes of the polymeric 
cable materials, which is beyond the scope of this paper. 

II. HEAT CONDUCTION MODEL OF A SINGLE CABLE 

, A cable consists of a core containing the insulated wire conductors 
and an outer protective sheath. In general, the cable core will contain 
a loose array of copper wires (which constitute 40 to 50 percent of the 
core) covered with a thin, polymeric insulating layer. The dominant 
path of heat conduction, because of their high conductivity, is in the 
longitudinal direction of the copper wires. Heat is also transferred 
radially, through the porous array of wires, by virtue of thermal 
radiation and heat conduction through contact points along the length 
of the wires, as shown in Fig. 2. The cable core cannot be considered 
a continuum because of the noncontiguous nature of the wires; hence, 
in developing the governing equation of heat conduction it is desirable 
to consider the individual wires as microstructural elements. This 
concept, as will be seen, permits interaction of the wires and leads, in 



* Standard Methods of Fire Tests of Building Construction and Materials, American 
Society of Testing Materials. 
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Fig. 2 — Heat transfer modes in firestopped cables. 

the limit, to an approximate continuum representation of heat con- 
duction in the cable core. This development follows. 

Consider the individual conductor wires as composite cylinders with 
copper wire radius r c (see the appendix for a list of parameters) on 
which is affixed an insulating material with thickness (r d - r c ), as 
shown in Fig. 3. When radial symmetry prevails, Fourier heat conduc- 
tion equations and boundary conditions for the conductor wire assume 
the following forms: 



ra 2 T c 
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l d ( dTcu\ 

r ir V dr ) 
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Fig. 3 — (a) Macrocoordinate and (b) microcoordinate systems. 

for the insulation, with boundary conditions: 
T cu (0, z, t) = finite, 
T cu (/y, z, t) = Ti(r c , z, t)* 

dT cu dTi 

k cu — — {r c , z, t) =ki-— (r c , z, t), 
Br dr 

BT- 

hi — (r d , z, t) = -hi[Ti{r dy z, t) - f(z, t)], 
Br 



and initial conditions: 



T cu (r, z, 0) = Tdr, z, 0) = 0. 



(3) 
(4) 

(5) 
(fi) 

(7) 



* A perfect contact between the copper wire and insulation is assumed, since the 
insulation is literally molded onto the wire. 
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In the above equations, T(r, z, t) is the temperature, k is the thermal 
conductivity, pc is the volumetric sensible heat capacity, z is the axial 
coordinate along the length of wire, r is the radial coordinate, t is the 
time, hi is an effective linear heat transfer coefficient, and T(z, t) is a 
temperature to be assigned subsequently. The subscripts "cu" and "i" 
are assigned to quantities associated with the copper conductor and 
insulation, respectively. Since the diameter of the wire is much smaller 
than its length, it is convenient to express the temperature in the series 

T cu (r, z, t) = T (z, t) + T™(r, z, t) + • • -, (8) 

and 

Ti(r, z t t) = T (z, t) + Tl l) (r, z, t) + ■ ■ ■ . 0) 

The leading term in eqs. (8) and (9) is the average temperature, i.e., 



T (z, t)= — 



f 

Jo 



rT cu (r, z, t)dr + 



rTi(r, z t 



t)dr 



(10) 



from which it follows that 






rT?J{r, z, t)dr + 



J 



rT} n) (r, z, t)dr = 0, 



n — 1, 2, 3, 



(11) 



Substituting series (8) and (9) into differential equations (1) and (2) 
produces a set of recurrent differential equations: 



Bz 2 r dr \ dr / 
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n = 1, 2, 



(12) 



(13) 



All TlVir, z, t) and T\ n) {r, z, t)n = 1, 2, 3, ... are then expressible in 
terms of the average temperature, T (z, t). Retaining terms up to first 
order, it follows from (12) and (13) that 

r 2 
T cu (r, 2, t) = T (z, t) ~j^L cu + A(z, £)ln r + B(z, t), (14) 



and 



where 



Ti(r, z, t) = T a (z, t)- — V + C(z, t)\n r + D(z, t), (15) 

4/2, 
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The variables of integration A, B, C, D are determined from boundary 
conditions, eqs. (3) through (5) and eq. (11). The final form of the 
temperature in the copper wire and insulation in terms of the average 
temperature T (z, t) is 

T cu (r, z, t) = T (z, t) 

+ i( r? -fi-' 2 ) L ~ + l ln( '^ )(LC "- L ' ) 



+ 



4kird 



-rlU 



+ W±** 



(16) 



and 



, r^ln(r/r d ) 
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8fe 






(17) 



It now remains to obtain the differential equation for the average 
temperature, T ((z, t). This equation comes from the satisfaction of 
the boundary condition at r = r</, as shown in eq. (6); 



1 r c fli 2 2\ 



L l + 



2r 2 c (r 2 d - rl) 



frt_hi_ 

\u u = y d [T {z,t)-f{z,t)l (18) 



Relating the variables in each conductor wire to those of neighboring 
ones in the limit leads to a continuum base from a discrete one. To 
arrive at this continuum description of the cable core consider the 
wires as embedded in a radially symmetric macrocoordinate system, 
as shown in Fig. 3. The center of, say, the nth conductor is located at 
R n and eq. (18) is written at this point. In keeping with the assumption 
of radial symmetry, the temperature on the right of (18), T(z, t), is 
related to the average temperature of adjacent wires T { " +1} and T£" _1) 
located at R n + 2r d and R n — 2/v, respectively. To determine T(z } t) 
the temperatures are weighted with respect to the location of the wire 
in the cable core. This has the effect of enforcing radial symmetry 



f(z, t)=-^r- [(ft. + r d )T ( r l) + (Rn - r d )T<"- 1) ]. 

Attn 



(19) 



Thus, if D is used to denote the difference appearing on the right of 
(18), it can be written with the aid of (19) as 
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D = -^~ [<*» + r d )T { r l) - 2RnT™ + (R n - r d )T { n ~ l) l (20) 

When we expand T { <? +1) and T^ _1> in a Taylor series about R n , eq. 
(20) becomes 

— ¥(i [*£]♦•»}■ » 

If we substitute (21) into (18) and consider r d <K 1 so that only the 
first term on the right of (21) is significant, the final form of the 
continuum representation of the cable core, after the subscript n is 
eliminated, is obtained as: 
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The quantity (4A,r d ) has the dimensions of conductivity and represents 
the effective radial conductivity of the cable core. The conductivity in 
the axial direction reflects the effect of the microstructure and is given 
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and the effective heat capacity as 
fii(r d — r c ) 
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Notice that when the conductor wires are assumed insulated, i.e., 
hi = 0, there is no radial heat flow and the effective heat capacity and 
axial conductivity are given, as expected, by the law of mixtures. 

Now that the heat conduction equation for the cable core has been 
determined, it remains to obtain the equation that governs the heat 
transfer in the entire cable including the outer sheath. This develop- 
ment parallels the procedure employed above. The only exception is 
that, since the cable core and sheath are not generally in perfect 
contact, the average temperatures in the cable core and sheath are 
assumed to be different. In many cables the cable core can be pushed 
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through the cable sheath with only a moderate pressure. Since the 
axial conductivity of the cable core is much larger than that of the 
cable sheath, as a simplification, it appears reasonable to assume that 
the axial heat conduction in the cable sheath can be neglected. Thus, 
the expression for the temperature in the cable core and sheath 
assumes the form 



T (R, z, t) = 7T(z, t) + {2R \ * g) L™(z, t), 



and 
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where 
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and Tf\z, t), T s Q) {z, t) are the average cable core and cable sheath 
temperatures, respectively. The other parameters are defined in the 
appendix. The average cable core and cable sheath temperatures given 
above are related through the imperfect heat transfer boundary con- 
dition at R = R c , namely, 
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(z, t) + h ca T™(z, t), (26) 



where h ca is the heat transfer coefficient between the cable core and 
cable sheath. 
In the following section eqs. (24b) and (26) are used to determine 
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the heat transfer in an array of sheathed cables that form a bundle. 
However, before this is accomplished it is worthwhile to briefly sum- 
marize the preceding developments. 

The cable core is obviously not a continuum; consequently, a micro- 
structural characterization starting with an individual conductor wire 
was used to derive the governing radially symmetric heat conduction 
equation for the cable core. This derivation produces an effective or 
averaged macroscopic continuum representation of the cable core that 
retains the physical characteristics of the conductor wire microstruc- 
ture. It is, however, anisotropic, since the thermal conductivity in the 
axial and radial directions is not equal. After we obtain the governing 
equations for heat flow in the core, the equations that the cable sheath 
temperature satisfies [eqs. (24b) and (26)] were then determined. 
These were obtained, paralleling the preceding microstructural deri- 
vation, for heat flow predominantly in the axial direction. The radial 
heat transfer is introduced as a perturbation. 

III. HEAT TRANSFER THROUGH A CABLE BUNDLE 

The objective of this investigation is to analyze heat flow through 
firestopped cable bundles during fire tests. The ASTM El 19 temper- 
ature-time history, Fig. 1, provides the fire environment. A typical 
firestopped cable bundle configuration and the above- and below-floor 
coordinate system is illustrated in Fig. 4. 

As previously mentioned, cable penetrations vary in size and in the 
number and type of cables accommodated; to lend a degree of defi- 
niteness to the analysis it is convenient to consider a widely used 
arrangement. A square cable array containing nine cables, as illustrated 
in Fig. 5, is a mathematically manageable configuration, yet most of 
the heat transfer characteristics of larger cable arrays are maintained. 
From the point of view of symmetry, only three cables need be 
considered — the corner cable, to be designated hereafter as cable 1, 
the side cable, to be denoted as cable 2, and the center cable, to be 
called cable 3. 

The section of each cable below the slab (see Fig. 4) is directly 
exposed in a furnace to a fire temperature of up to 1000 degrees 
centigrade and, therefore, after only a short time all polymeric insu- 
lating material is burned away. In the model this effect is approximated 
by assuming that only a loose array of independent copper wires 
projects below the slab. This section of the cable bundle extending up 
through and above the floor slab experiences a different thermal 
environment. A zone of decomposition of the polymetric insulating 
materials occurs and creeps upwards during the extended period of 
exposure to the below-floor fire. This zone of charred and expanded 
insulating material alters the temperature in the void spaces. In the 
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Fig. 4 — Typical cable opening firestop test configuration. 



present analysis the effect of this nonlinear phenomenon is ignored to 
render the analysis tractable. It is conservatively assumed here that 
the hot furnace gases move unimpeded in the void space between 
cables. 

The average cable core and average cable sheath temperature for 
each of three cables is determined by performing a heat balance at the 
outer cable sheath surface R — R s . This leads to the following three 
coupled equations: 



dT 



a) 



(0) 



dR 



(R a , z, t) = 2A c Hcb[T™ - Ti l} ] + (3A f + 2A c )H f [Tr - T™] 
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Fig. 5— Cable bundle configuration considered. 
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dR 



+ AfHtiT? - Ti 1 *] + 2A f HrF 12 [_T™ - T' 1 '] 

+ A f H r F 2 2[Tf ) - T< u ] | r-r, (27a) 

(R S} z, t) = 2Atf cB \T™ - Tf] + (AMcb + 2A f H r F^ 



. [T ( 3 ) _ T ( 2)] + 2A f H 1 [Tf ) - T™] + (2A f +A c )H f [T ( ; ) - T? ] 

+ 2A,H P F ia [Ti lt -T. (2) ]| M . (27b) 

k B ^-(R s , z, t) = AA c H cB [T ( P - T« + 4A,/f 1 [T} 1) - Tf ] 

Oil 

+ S^H^fTf - T< 3) ] + 4A f H r F 22 [Ti 1) - T™]\ r =r,. (27c) 

In the above equations the parameter A c denotes the fraction of the 
cable sheath perimeter that is in contact with an adjacent cable, and 
A { the portion exposed to the hot furnace gases that flow up through 
the interstices of the cable. It therefore follows that A{A C + A f ) = 1. 
H cB is the solid-contact conductance between the cables, Hf is the heat 
transfer coefficient between the cables and the firestop material, Hi is 
the convective heat transfer coefficient between the cable and the hot 
furnace gases that flow up through the void spaces between cables, 
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and H r is a linearized black body radiative heat transfer coefficient 
between the cables. F i2 and F 22 are radiation form factors between 
cable 3 and cable 2, and cable 1 and cable 3, respectively. The 
temperatures TfHz, t) and T l f 1} (z, t) are assumed known and represent 
the temperature distribution in the mineral wool insulating flrestop 
material surrounding the cable bundle and the temperature of the 
furnace gases that move up through the void spaces between cables, 
respectively. 

At this juncture, it is convenient to eliminate the independent time 
variable t from the equations by introducing the Laplace transforma- 
tion 



T(z;p)=\ T(z, t)e~ pt dt. (28) 



>AP)= \ 
Jo 



In what follows, a bar over a variable indicates that the transformation 
(28) has been performed. Substituting eq. (24b) into (27) and using eq. 
(26) to eliminate the average cable sheath temperature, we produce a 
system of three coupled ordinary differential equations for the trans- 
formed average cable core temperatures f o?i , ?$, and T$l ' 

(r,£> 2 - Ax)f a - 2( yi D 2 - ft)!T8 - {ijiD 2 - eOTB 
= -(ha, + pEJUSAf + 2A c )H f ff ) + A f H,f { Pl (29a) 

-2( Y i£ 2 -0i)f o } + (V 2 D 2 - A 2 )T$ - (r, 2 D 2 - e 2 )f$ 
= -(h ct + pE 1 )[(2A r + A c )H f f f + 2A f H 1 f i f l) l (29b) 

-M-mD 2 - ei)f 8 - (y^ 2 - &)?$ + (r 3 z> 2 - A*)T® 
= -4(h cs + pEMfHiff, (29c) 

d 2 
where the operator D 2 = -5-* and the coefficients T, A, y, 6, ij, and e are 

dz 

linear functions of the transform parameter p. After considerable 

algebra the system of differential equations (29) is solved for the three 

average cable core temperatures in the form 

Tf> = Ciexp(— z4m\) + C 2 exp(-zVm^) + C 3 exp(— z^/rm) 

-(h cs +pE l )\_RTf ) + HTfl (30) 

The quantities mi, m 2 , and m 3 are functions of the Laplace transform 
parameter p and are evaluated as the roots of the characteristic cubic 
equation of the system of equations in (29) . The coefficients of inte- 
gration Ci, C 2 , and C 3 are determined by matching the solution of the 
cable section below the slab. This is accomplished by enforcing tem- 
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perature and energy continuity for each of the three cables at the 
common boundary Z = and X = L, as shown in Fig. 4. The cable 
sheath temperature evaluated at R = R s follows from (30) in the form 



T a {R*,z,p) = 



h cs + pBi . . n 

- — ■ — =-(1 +pAi) +pDi 

flcs -r P&l 



1 



- (**), 



Bi(l+pAi) 



(hcs + pEi) 



+ D, 



d 2 T m 

^rf. (3D 

dz l 



where A\, J3i, D\ and E\ are constants. 

The temperatures presented in this section are in the Laplace 
transform domain and must be inverted to the real-time regime. The 
functions, however, are too complicated to be inverted in closed form. 
A numerical procedure using the method of quadratures to obtain 
these inversions is discussed in Section IV. 

IV. INVERSION OF LAPLACE TRANSFORM TEMPERATURE SOLUTION 

The form of the transformed temperature solution given in eqs. (30) 
and (31) is much too complicated to use for obtaining a closed-form 
inversion formula. Consequently, we must resort to a numerical inver- 
sion procedure. Most of the methods that appear in the literature 
involve expanding the transformed function in a series that could then 
be inverted term by term using tabulated formulae. Littlewood and 
Zakian 1 suggest expanding in a series of Chebyshev polynomials, while 
Longman 2 proposes using the Pade table for the Taylor series expan- 
sion of the transformed function. Both of these methods were judged 
to be impractical because of the complexity of the functions to be 
inverted here. 

The method that was finally adopted was developed by Talbot. 3 The 
inversion of arbitrary transforms is accomplished by a method of 
quadrature along a special contour in the complex plane. The standard 
inversion formula for a transformed functionF(p) involves performing 
the following integration in the complex plane: 

F(t)=—. F(p)e pt dp, (32) 

2m J . 

where p is considered a complex variable, i = (-1) 1/2 , and y is to the 
right of all the singularities ofF(p). In Talbot's method the path of 
integration indicated by (32) is deformed to the path L shown in Fig. 
6, the equation of which is 

p = a + i$, (33) 
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Fig. 6 — Path L for Laplace transform inversion. 

where a — 6 cot 6. 

Path L is equivalent to the standard integration path providing that 
(i) L encloses all singularities of F(p), and 

(ii) | F(p) | — ► uniformly mRep<0as\p\ — » oo. 
Condition («) holds for the functions considered here. Condition (i), 
in general, may not be satisfied by a given F(p); however, the modified 
function F(Xp + a) for suitable scaling and shift factors \ and a can be 
made to conform. In this regard, for the functions to be inverted here 
it can be argued that the singularities exhibited by eqs. (30) and (31) 
are located on the negative real axis; consequently, condition (i) is 
satisfied without resorting to scaling or shift factors. However, to 
accurately perform the inversions over large intervals of time (the 
temperatures are calculated for times up to two hours) a scaling factor 
A > 1 is necessary. This scaling factor merely shifts the singularity 
along the negative real axis closer to the origin. 

Once conditions (i) and (ii) are satisfied, the inversion formula (32), 
when taken about L assumes the form 



Fit) 



2iti J 



e" T F(\p + o)^d8, 



r = Xt 



(34) 



Finally, "M>" Simpson quadrature with equal intervals w/n in the 
variable gives the approximation 
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Ap ot 

F( t ) {Re{H + 4Hi + 2H 2 + 4H 3 + ■ ■ ■ + 4tf„-i)}, (35) 

3n 

where 

ff*«ff(p)U, 

Ok = Aw/" <^ = °> •••»»— 1) 
H(p) - e PT F(Ap + a)(l + #), 
and 

- [fl - % sin(20)]/sin 2 0. 

The symbol ite indicates that only the real part of the complex 
quantity is taken. Sufficient accuracy is obtained by a suitable choice 
of n t X, and a. The principles of choice are presented by Talbot 3 and 
will not be discussed here. 

A check of this technique was made by inverting the function given 
in Ref. 2. In all instances satisfactory results were obtained. For 
numerical inversions performed here, it was found sufficient to take 
n = 20, \ = 8, and a = 0. 

V. DISCUSSION AND SUMMARY 

The concept of firestopping any type of penetration is synonymous 
with retarding the flow of heat from the fire side of the penetration to 
the unexposed side. The problem rests solely on the identification of 
procedures and materials to seal the space adjacent to the penetrant 
to meet certain standards. Cables by their very nature have excellent 
thermal conduction properties in the longitudinal direction. Very little 
can be done to prevent heat conduction up the array of copper wire 
conductors that make up the center core section of the cable. However, 
it appears reasonable to measure the efficiency of a firestop by the 
temperature rise on the free surfaces on the unexposed side of the 
penetration. If these surfaces are maintained at sufficiently low tem- 
peratures, ignition of combustibles that happen to be in direct contact 
or in close proximity cannot occur. In the case of cable penetrations 
the critical surfaces are the horizontal firestopping material surface on 
the unexposed side and the vertical cable sheath surfaces at the outer 
perimeter of the cable bundle, as shown in Fig. 4. In normal practice 
the firestop material is chosen to be a good thermal insulator and is 
applied with sufficient depth to preclude high temperatures at the top 
surface. Thus, the efficiency of the firestopped geometry will, in 
general, be determined more by the temperature rise on the cable 
sheath surface of the outer perimeter cables than by the firestop 
material. To ascertain the temperature at this critical location, a 

TRANSIENT HEAT TRANSFER 2327 



theoretical transient heat transfer model of a firestopped cable bundle 
was developed. In addition to heat flow through the cable core, the 
model treats the lateral heat flow from the cable core to the cable 
sheath and convective heat transfer to the cable sheath from the flow 
of combustion gases through the interstices between the cables. The 
standard ASTM El 19 temperature variation is applied at one end of 
the cable bundle and the temperature distribution of the individual 
cables is computed along the length for up to two hours. To quantify 
these effects, temperature distribution was computed at 10-minute 
intervals for two hours at the 20-cm and 30-cm firestopped depth for 
three different cable types — switchboard, terminating, and power ca- 
ble. For each of these cable types the following parameters were 
varied: 

(i) R s , cable radius 

(ii) Ac, fraction of cable sheath surface in contact with an adjacent 
cable 

(iii) h cs , conductance between cable core and sheath 

(iv) Hf r conductance between cable bundle perimeter and firestop 
material 

(v) Ap, furnace pressure. 
Some remarks concerning the inclusion of the furnace pressure as a 
parameter are in order. Furnace pressure influences the magnitude of 
the heat transmitted to the cable sheath — by the hot furnace gases 
that travel up through the space between cables — through the heat 
transfer coefficient Hi, as shown in eq. (27). The value of Hi is 
computed from standard empirical formulae 4 once the gas flow velocity 
and flow-channel characteristics are known. The steady-state gas 
velocity is calculated from the furnace pressure assuming that the 
spaces between cables are independent flow channels. 5 In general, Hi 
increases with increased furnace gas pressure. 

Some of the results of the analysis for a given type cable are given 
in Table I and Fig. 7. The calculated end-point temperature at z = 30 
cm in Fig. 6 and t = 2 hours for the cable sheath and core are given in 
Table I, and a typical temporal temperature distribution is given in 
Fig. 7. The temperatures given in Table I and Fig. 7 should not be 
construed to be indicative of actual measured test values. As previously 
indicated, the model does not take into account nonlinear aspects of 
this obviously complex phenomenon, such as potential combustion 
modes and melting of the polymeric materials that in some instances 
could conceivably, for periods of time, constrict the void of spaces 
between cables and thereby reduce the flow of furnace gases. Never- 
theless, the sensitivity of the cable temperatures to changes in the 
linear cable bundle parameters identified in the model, H f , A C ,R S , and 
A/> can be calculated and are presented in Table I. 
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Table I — Calculated temperatures at z = 30 cm and f = 2 hours for 
various values of the cable bundle parameters 



(W/cm 2 - 
Case °C) 



A 



(cm) 



Cable Sheath 

Temperature 

(°C) Rise 



Cable Core 

Temperature 

(°C) Rise 



Ap (cm 
of Water) Cable 1 Cable 3 Cable 1 Cable 3 



0.034 
0.017 
0.017 
0.017 

0.017 



0.1 
0.1 
0.15 
0.1 

0.1 



1.27 

1.27 

1.27 

0.635 

1.27 



1 
I 
1 
1 
0.13 
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224 
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93 
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Fig. 7— Cable core temperature for Case 1 cable bundle parameters given in Table I. 

It is clearly shown in Table I that the cable core temperature exceeds 
that of the cable sheath. A tightly packed firestop that exerts lateral 
pressure on the side of the cable array will provide sufficient heat 
sinking to reduce the temperature in the cable bundle. This physical 
effect is embodied in the contact conductance parameter, Hf. The 
larger the numerical value of this parameter, the tighter the firestop 
packing. Cases 1 and 2 of Table I show that increasing this parameter 
indeed results in a lowering of the cable sheath and cable core tem- 
peratures. 

The same general result prevails when furnace overpressure An is 
reduced. This effectively reduces the heat transfer coefficient Hi 
between the cable sheath and the hot furnace gases. This is demon- 
strated in Table I by comparing the end-point cable temperatures of 
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Cases 2 and 5. Cable temperature can also be reduced by increasing 
the contact surface between cables, as measured by an increase in the 
parameter A c , suggesting a tighter and more compact cable bundle. 
The cable temperatures, as observed in Cases 2 and 3, are also lowered. 
This is primarily due to the resulting smaller void space between the 
cables and secondarily to the larger conduction path presented to the 
interior cables. 

VI. CONCLUSIONS 

The following trends were generally indicated: 
(i) The cable core and cable sheath temperature is largest for the 
interior cable (cable 3 of Fig. 5) and the smallest for the corner cable 
(cable 1 of Fig. 5). 

(«) The temperature of cable core exceeds that of the cable sheath. 

(Hi) The primary heat transfer mode to the cable sheath is from 
the flow of hot combustion gases through the void space between 
cables. 

(iv) Reducing the void space between the cables by tightly packing 
the cables and/or using smaller diameter cables impedes the flow of 
hot combustion gases and results in a significant reduction of the 
primary convective heat transmission mode. 

(v) A tightly packed firestop capable of providing some heat sink- 
ing to the cooler environs is the most practical and effective method of 
reducing the heat transfer properties of the cable bundle. 

(vi) The magnitude of the cable sheath temperature for similar size 
cables depends on the furnace gas pressure and, to a lesser degree, on 
the firestop depth. 
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APPENDIX 

A c — Fraction of cable sheath surface in contact 

with an adjacent cable. 
A f — Fraction of cable surface in contact with com- 
bustion gases, 4(A? + Af) = 1. 
jf l2i F 22 — Radiation form factors between cables. 

" Hi— Heat transfer coefficient between cable and 

combustion gases. 
Hcb— Solid contact conductance between cables. 
/^—Contact conductance between cable bundle 

perimeter and firestop. 
H r — Linearized black body radiative heat transfer 

coefficient. 
hcs— Conductance between cable core and cable 

sheath. 
A,— Conductance between individual copper con- 
ductor wires. 
k£ = AhiTd— Effective radial thermal conductivity of cable 
core. 
k eu — Thermal conductivity of copper wire. 
ki— Thermal conductivity of wire insulation. 
jfe s — Thermal conductivity of cable sheath. 
(A z )eff— Effective axial thermal conductivity of cable 
core. 
Kcu = k cu j{pc) cu — Thermal diffusivity of copper wire. 
Ki = ki/ (pc)i— Thermal diffusivity of wire insulation. 
k s = k s /{pc) x — Thermal diffusivity of cable sheath. 
(K g U={k,U/(pcU— Effective axial thermal diffusivity of cable 
core. 
L— Length of cable below slab. 
R c — Radius of cable core. 
# s — Overall radius of cable. 
r c — Radius of single copper conductor wire. 
r d — Total radius of single conductor wire. 
Ap — Furnace pressure. 
(pc) C u— Sensible volumetric heat capacity of copper 

wire. 
(pc) e ff— Effective sensible volumetric heat capacity of 
cable core. 
(pc)i— Sensible volumetric heat capacity of wire in- 
sulation. 
( p c) s — Sensible volumetric heat capacity of cable 



sheath. 
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